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ABSTRACT 
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.'This  paper  examines  the  concepts  of  feedback  and  Adaptivity 
for  the  Finite  Element  Method.  The  model  problem  concerns  P  ' 
elements  of  arbitrary,  fixed  degree  for  a  one-dimensional  two-point 
boundary  value  problem.  Three  different  feedback  methods  are  intro¬ 
duced  and  a  detailed  analysis  of  their  adaptivity  is  given. 


1.  INTRODUCTION 


In  recent  years  there  has  been  a  growing  interest  in  designing 
"adaptive"  methods  for  the  numerical  solution  of  differential  equations. 
"Adaptive"  methods  for  numerical  quadrature  and  integration  of  initial 
value  problems  for  ODE's  are  common  in  available  codes.  The  term  "adaptive" 
however  is  generally  used  in  a  very  loose  sense. 

In  this  paper  we  shall  attempt  to  give  a  more  precise  analysis  of  the 
concepts  of  feedback  and  adaptivity  as  they  relate  to  the  numerical  solu¬ 
tion  of  differential  equations.  The  definitions  of  feedback  and  adaptivity 
that  we  use  are  similar  to  those  found  in  [11]  (see  also  [  7  ]):  a  feedback 
method  is  characterized  by  the  fact  that  it  produces  a  sequence  (a  tra¬ 
jectory)  of  approximate  solutions,  with  each  new  solution  depending  on 
information  available  from  the  preceding  solutions;  a  feedback  method  is 
called  adaptive  if  it  is  optimal  with  respect  to  some  performance  measure, 
more  precisely  if  its  performance  with  respect  to  this  measure  is  not 
worse  than  that  of  any  other  trajectory. 

We  restrict  the  discussion  to  finite  element  methods  for  a  one¬ 
dimensional  two-point  boundary  value  problem;  the  elements  are  C°  piece- 
wise  polynomials  of  arbitrary,  fixed  degree.  The  performance  measures  we 
consider  are  related  to  the  convergence,  respectively  the  quasi-optimality, 
of  the  approximate  solutions.  We  believe  the  approaches  and  the  theorems 
discussed  in  this  paper  to  be  valid  in  dimensions  higher  than  one;  a  feed¬ 
back  method  of  the  same  type  as  those  introduced  here  has  been  implemented 
for  a  second  order  elliptic  system  of  two  equations  on  a  two  dimensional 
domain  in  the  program  FEARS  (cf.  [10]).  Both  the  available  theoretical 
results  (cf.  [2]  [  3  ]  [  6  ]  and  the  extensive  numerical  experimentation 
strongly  support  our  belief. 
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The  feedback  methods  constructed  here  show  some  similarity  to  the 
algorithm  for  approximation  introduced  in  [12],  as  well  as  the  ODE-solver 
studied  in  [  9  ]  in  the  sense  that  they  try  to  achieve  an  equilibration  of 
certain  (locally  computable)  quantities  associated  with  the  error.  See 
also  [  8]  [13]  for  a  treatment  of  elliptic  P.D.E. 

The  organization  of  this  paper  is  as  follows:  In  section  2  we  intro¬ 
duce  the  model  two  point  boundary  value  problem,  its  finite  element  solu¬ 
tion  and  define  the  error- indicators  and  the  error  estimator  on  which  all 
the  feedback  methods  are  based;  the  definitions  are  similar  to  those  in 
[  5  ] •  Sections  3  through  5  make  precise  our  notions  of  feedback  and  adap¬ 
tivity  and  provide  three  examples  of  practically  interesting  feedback 
finite  element  methods.  Theorem  6.1  is  our  first  result  concerning  adap¬ 
tivity  of  feedback  methods,  it  shows  that  for  a  very  large  class  of  feed¬ 
back  methods  the  corresponding  approximate  solutions  will  always  converge 
in  energy  to  the  exact  solution.  In  section  7  we  examine  the  question  of 
adaptivity  with  respect  to  the  stronger  performance  measure  of  quasi¬ 
optimality.  The  Theorems  7.1  and  7.2  verify  that  the  three  examples  from 
sections  4  and  5  all  produce  approximate  solutions  that  are  quasi-optimal , 
provided  the  exact  solution  satisfies  some  additional  assumptions;  we  go 
to  some  length  to  show  that  these  assumptions  are  very  mild  and  "almost" 
necessary.  In  order  to  compare  two  different  feedback  methods  that  are 
both  adaptive,  a  relevant  quantity  is  the  cost  of  the  computation.  In 
section  8  we  give  a  very  simple-minded  comparison  of  the  feedback  methods 
introduced  in  Sections  4  and  5.  Finally  the  appendix  contains  the  proof 
of  a  lemma  that  establishes  the  equivalence  of  our  error  indicators  and 
a  local  interpolation  error.  This  lemma  was  used  in  verifying  quasi-optimality 
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and  the  proof  of  it  is  based  on  a  stability  estimate.  Lemma  9. 
is  a  slight  extension  of  a  result  contained  in  [4]. 


1,  which 


4 


2 .  PRELIMINARIES 

Consider  the  two-point  boundary  value  problem 


(2.1a) 

Lu  = 

~h a(x) 

du 

dx 

+  b(x) 

u  =  i 

:  in  7 

(2.1b) 

u(0)  = 

u(l)  = 

0, 

» 

where 

a,b  C  L  (7) 

00 

with 

(2.2a) 

0 

al  S 

a  (x) 

~  c2 

(2.2b) 

0 

< 

b(x) 

'  a2  ’ 

and 

(2.3) 

f  €  H~ 

X(7). 

The  solution  of  (2.1a-b)  is  to  be  understood  in  the  usual  weak  form, 
i.e.,  u  €  H^d)  and 


B(u,v) 


f1 

^0 


fvdx 


for  any  v  €  Hu),  where  B  denotes  the  bilinear  form 


(2.4) 


B(u,v)  = 


(a  —■  4^  +  buv)dx. 
dx  dx 


°1 

The  space  H  (T)  consists  of  all  functions  that  have  first  derivates  in 

-1  °1 

L2O )  and  that  vanish  at  0  and  1.  H  (7)  denotes  the  dual  of  H  (I). 
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We  define  the  energy  norm 


Hull, 


=  !  B(u,u)  ] 


1/2 


For  any  mesh 


0  =  xl 


*N(A) 


=  1 


let 


hA  =  h(IA) 
3  3 


A  A 

x.  -  x.  , , 
3  J-l 


j  =  1,  .  .  .  ,N(£.) 


and  hu  =  max  h.;  x7  are  called  the  nodal  points,  it,  the  ele- 

j  2  2  2 

ments  and  N(j)  the  cardinality  of  the  mesh  A.  A  mesh  may  be  identified 

by  either  the  set  of  nodal  points  or  the  set  of  elements,  we  shall  thus 

write  A  =  {xA},  x  (  A  or  A  =  {1^},  I  (  A  depending  on  the  context. 

Let  p  2  1  be  an  integer.  By  S^(A)  we  denote  the  linear  subspace 
°1 

of  H  (I)  which  consists  of  functions  that  are  piecewise  polynomials  of 
degree  5  p  on  each  I  €  A. 

The  finite  element  solution  u(A)  t  Sp (A)  is  defined  in  the  usual 

way , 


(2.5) 


u(  A) 


P(A)u  , 


where  u  is  the  solution  to  (2.1a-b)  and  P(A)  is  the  elliptic  projec¬ 
tion  of  H^(I)  onto  S*\a)  .  By  this  we  mean  P(A)u  (  S^(A)  and 


(2.6) 


B(P(A)u,v) 


B(u,v) 
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p 

for  any  v  (  S  (A) .  The  expression 

(2.7)  e ( A)  =  u  -  u(A) 
is  the  error  of  the  finite  element  solution. 

ol 

Let  R(A)  and  R (I \I)  be  the  elliptic  projections  of  H  (I )  onto  the 
spaces 

{v  t  H"^ (7  )  ;  v(xT')  =  0,  j  =  1, . .  .  ,N(A)  } 

and 

{v  f  H1(I)  :  v (x)  =  0,  x  €  I  \l} 

respectively,  with  this  notation 

N(A) 

(2.8)  R  ( A )  v  =  l  R(l\l?)v. 

j-1  J 

We  define  the  error  indicators  y^(l7)  =  and  the  error  estimator 


t(A)  as  follows 

(2.9) 

A 

"j  = 

|lR(J\Ij  )e(A)  1 

(2.10) 

f(A)  = 

1 R(A)e(A)  |'E  . 

Due  to  (2.8)  we  conclude  that 


(2.11) 


The  functions 


V4) 


RO\I?)e(M, 


which  enter  into  the  definition 


(2.9)  of  the  error  indicators,  are  characterized  by 


(2.12a) 

Lw.(A)  =  rl  =  f  -  Lu(i) 

J  3 

in  IT 
3 

(2.12b) 

w  (a)  =  0  in  I  \  1 1  . 

To  simplify  the  analysis  in  this  paper  we  shall  assume  that  the  (£)'s 
are  directly  computable.  In  practice  one  derives  explicit,  approximate 

A 

formulae  based  on  the  residuals  rl  (cf.  [5]). 

J 

An  important  result  is  (see  [1]  [5]): 

Lemma  2.1.  There  exists  C  >  1  depending  on  and  a ^  of  (2.2a-b)  - 

but  independent  of  a,  b,  f  and  A  such  that 

(2.13)  t  (fi)  5  iie(A)i  F  S  Cf(A). 

For  completeness  we  provide  a  short  proof.  Let 


aA 

3 


a(x)dx. 


j  =  1,...,N(£) 


and  define 

(2.14) 


■SA(u,v)  = 


N(£) 

I 

3=1 


A 

Vi 


A  du  dv  , 
a  —  —  dx 
j  dx  dx 


on  H^(7)  v  H^(7).  P(A)  denotes  the  orthogonal  projection  of  H ^ ( 7 )  onto 
SP(a)  using  the  inner-product  (2.14).  It  is  well  known  that 


P(j,)v(x“)  =  v(xi'). 


for  any 


v  i  H1(I). 


Thus 


j  =  0,  .  .  .  ,K(:.) 


(2.15) 


e(A)  E  =  B(e(fl),e(A>) 

=  B(e(A),e(i)  -  P(A)e(/») 

=  B(R(A)e(A),e(A)  -  P(A)e(A)) 

5  C,'  R(A)e(h)  i  £  e(l)  1  E  . 


In  (2.15)  we  used  the  fact  that 

e(A)  -  P(A)e(l)  =  0  at  the  nodal  points, 

we  furtheremore  used  the  equivalence  of  the  energy  norms  associated  to  B 
and  Ba: 

£  B,(v,v)  <  B(v,v)  5  CB  A  (v ,v) 

with  C  only  depending  on  and  c  .  The  inequality  (2.15)  immediately 

leads  to  the  last  inequality  of  (2.13).  The  first  inequality  in  (2.13) 
follows  from  the  fact  that  R(£)  is  an  orthogonal  projection  in  the  inner- 
product  B.  a 


With  additional  assumptions  on  the  coefficients  a  and  b  we  may 
prove  that  not  only  is  the  estimator  equivalent  to  the  finite  element  error 
it  is  indeed  asymptotically  exactly  equal  to  the  error. 

Lemma  2.2.  If  in  addition  to  satisfying  (2.2a-b)  a  and  b  have  the 


i  f 


property  that 


9 


(d^}  3’ 


(j-)k_1  b 
dx 


L  (7) 


for  some  integer  k  >  1,  then 

e(i).£  =  r(A)(l  +  0(h")2~) 

'  t 

with  (j.  =  min(p,k)  .  The  term  0(h")  is  bounded  by  Ch“  independently 
of  f  and  L . 

Proof.  Let  R^Ci)  be  the  elliptic  projection  onto  the  orthogonal  com¬ 
plement  of  R(A)H^(7).  Then 


e(d)  =  R(d)e(£)  +  RJ(d)e(') 


^(A)  +  e?(i) , 


and 


(2.16)  Le^CA)  =  0  on  each  element  if. 

From  (2.16)  and  the  regularity  assumptions  about  a  and  b  it  follows 
that 


e2  I  k+1  > 

2  Hk+1(lJ) 


c'e2(  A)  x 

H  (Ip 


where  H  denotes  the  standard  Sobolev  space  of  functions  with  derivatives 
up  to  and  including  order  k  in  L2(lf)  (the  constant  C  is  independent 
of  if).  Thus  there  exists  v  €  SP(fi)  such  that 


62(A)  =  v  +  z  and 


i!  z  1  . 

H 


(Ip 


C(hpH  e_(A)  l;  . 

1  2  H^l6) 

J 


with  p  =  min(p ,k) , 
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in  particular  we  see  that 

(2.17)  z  <  C(h^)H:e_(£)'  . 

E  2  E 

Since  0  =  B(e(A),v)  =  B(e^(£),v)  +  B(e2(A),  we  compute  that 
B(e2(i) ,e2(A) )  =  B(e2(A),v)  +B(e2(£),z) 

=  -B (e^ ( A)  ,  v)  +  B(e2(A),z) 

=  -B (e^ (a)  , &2 (Z.) )  +  B(e1(A),z) 

+  B(e2(£),z) 

=  Bte^A)  ,z)  +  B(e2(A)  ,z) 

(in  the  last  equality  we  used  that  e^(A)  and  e^(L)  are  orthogonal). 
As  a  consequence  of  this  and  (2.17)  we  obtain 

!-e2(A)l!j  £  Ci  e;L(A)i  E  e2(A)  E(h^)“  +  C  e2(A)  ^(h:)“, 

and  so  for  h^  sufficiently  small 

!!e2(A)!iE  £  C|!e1(A)'E(hi)u  . 

This  leads  to 

!!e(A)||E  =  !|e1(A)t;^  +  e2(£):^ 

£  ||ei(A)i!E(l  +  O(hV^) 

-  e(A)(1  +  O(hV^), 

thus  proving  the  desired  result  for  sufficiently  small.  For  large 

A 

h  this  lemma  is  a  simple  consequence  of  Lemma  2.1.  c 
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3.  FEEDBACK  AND  ADAPTIVITY  FOR  FINITE  ELEMENT  METHODS 

A  standard  finite  element  method  is  based  on  a  single  given  mesh  , 
whereas  a  feedback  finite  element  method  (feedback  f.e.m.)  produces  a 
sequence  of  meshes  A  and  corresponding  finite  element  solutions  u^, 
i  =  1,2,...  .  The  construction  of  the  meshes  L  is  governed  by  a 
transition  operator  A. 


"i+1 


A(A, 


'■i;u]L, . . .  ,ui) 


A  feedback  f.e.m.  may  be  terminated  through  some  stopping  criterion,  e.g., 
A  ,  u^  are  computed  until 

(3.1)  reap  5  -“Ve 


where  -  is  a  given  tolerance  and  c(a)  is  the  error  estimator  from 
before . 

If  a  measure  of  performance  is  provided  then  a  feedback  f.e.m.  shall 
be  called  adaptive  if  it  is  optimal  with  respect  to  this  measure.  We  do 
not  elaborate  on  the  most  general  definitions  of  feedback,  performance 
measure  and  adaptivity,  instead  we  refer  the  interested  reader  to  f 1 l ] . 

In  this  paper  we  analyze  some  concrete  methods  related  to  the  problem 
(2.1a-b)  and  our  performance  measures  evaluate  the  asymptotic  conver¬ 
gence  of  the  u^ .  In  practice,  when  analyzing  the  performance  of  a  feed¬ 
back  f.e.m.,  a  most  relevant  measure  is  the  computational  cost  in  satisfy¬ 
ing  a  certain  stopping  criterion — this  item  is  briefly  discussed  in  sec¬ 
tion  8.  One  important  question  that  enters  into  the  evaluation  of  the 
performance  of  a  terminated  feedback  f.e.m.  is  the  effectivity  of  the 
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stopping  criterion.  We  shall  not  address  this  issue  in  the  present  paper, 
but  only  point  out  that  the  result  of  Lemma  2.2  is  very  crucial  to  the 
effectivity  of  the  stopping  criterion  (3.1),  based  on  r ( A ) . 

For  simplicity  we  consider  only  binary  meshes,  i  .e. ,  the  nodal  points 
all  have  the  form  k/2J  for  some  0  5  j,  0  5  k  5  2J .  The  transition 
operators  never  delete  nodal  points,  so  that  the  meshes  L ^  form  an  in¬ 
creasing  sequence.  We  shall  call  a  transition  operator  simple  if 
only  depends  on  and  u^,  all  other  transition  operators  are  referred 

to  as  composite. 


"21 _ 
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4.  TWO  SIMPLE  FEEDBACK  FINITE  ELEMENT  METHODS 

t 

Let  and  let  nj^’  3  p  1,...,N(A^)  be  the  error  indica¬ 

tors  (2.9).  Below  are  the  definitions  of  two  very  useful  transition  opera¬ 
tors  and  A^;  in  each  case  we  list  the  sets  I  of  elements 

of  that  are  bisected  in  passing  to  a. 


..(1) 

“i 


1 


!_ 

S  max  n  (I)), 
If  A . 

l 


v<2) 

“i 


I(1) 

“i 


A.  , 

{I  1  }k 

1  j(*r  *«1 


where 


0  s  £  £  1. 


A,  A.  A. 

i  l  ^  i 

'j(l)  "  j(2)  "  -  nj(N(A1)) 

is  an  ordering  of  the  indicators  by 
size,  and  k  =  h*N(A.)] 

0  S  y  <  1.  1 

The  transition  operator  A^  bisects  all  elements  whose  indicators 

(2) 

are  larger  than  or  equal  to  S  times  the  largest  indicator.  A  bisects 
at  least  the  same  set  but  in  addition  always  bisects  the  [yN(A^)] 
elements  with  the  largest  indicators  (approximately  the  fraction  y) .  If 
3=0  both  of  these  transition  operators  bisect  all  elements.  Note  also 
that  A^  is  just  a  special  case  of  A^\  namely  corresponding  to 
y  =  0. 

There  is  of  course  no  reason  why  we  could  not  have  considered  dividing 
each  interval  into  n  *  2m  subintervals,  we  have  chosen  bisection  for 


simplicity . 
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5.  A  COMPOSITE  METHOD 

A .  1 . 

Let  I.  be  an  element  of  the  mesh  A..  A  direct  predecessor  of  I. 

2  1  A  2 

i 

is  an  element  I  £  A^,  1  s  k  <  i  -  1  with  the  property  that  1^  arises  as 
one  of  the  (two)  subintervals  in  a  bisection  of  I.  We  define 


A .  A. 

-  l  .  .  T  l 

I.  if  I.  €  Ax 


A  . 
1 


the  direct  predecessor  of  1^.  ,  otherwise. 


A.  Ai 

For  any  I.  £  A^  the  predictor  n(I^  )  is  given  by 

(A  A  .  „ 

[n  1(1.13  A.  A, 

— ^ — 3 — »  x-n  (i..  )j , 
n  k(I.) 

where  k  is  the  largest  integer,  Is  k  s  i,  such  that  I,  f  4,  and 

J  K 

0  <  x  1  is  a  constant.  The  f)  is  called  a  predictor  since  it  uses  past 

Ai 

experience  to  predict  the  effect  of  a  bisection  of  I_.  on  the  value  of 
the  indicator. 

(3) 

The  transition  operator  A  bisects  all  elements  of 


I^3)  :  {I.1  :  n  *(I  *)  »  6  max  fi(I)  > 

2  2  If A . 

i 

(3) 

for  some  fixed  OS?  i  1.  A  is  clearly  composite.  The  transition 
operator  used  in  the  FEARS  program  (cf.  [10])  is  a  direct  extension  of 
this  operator,  including  the  so-called  short  passes:  whenever  the  increase 
in  the  number  of  elements  is  less  than  a  preset  fraction  of  the  total 
number  of  elements  of  A^  no  finite  element  solution  is  computed  on 


are  constructed 


A.,,  .  Instead  indicators  for  the  elements  of  A.,, 
i+1  ,  l+l 

Ai  (3) 

directly  from  the  fl(I.  )  and  the  operator  A'  is  applied  to 

j  A  A 

A  A,, A  with  this  set  of  indicators  for  n  (I.  ).  This  process 

-L  *  •  •  •  9  l  1'tI  3 

continues  until  the  increase  in  the  number  of  elements  is  larger  than  the 


preset  fraction  of  the  total  number  of  elements  in  at  which  point  the 

short  pass  in  concluded  and  the  finite  element  solution  is  computed  on  the 
resulting  mesh  A^^  .  This  approach  avoids  the  computation  of  u(A^)  after 
only  a  small  increase  in  the  number  of  elements  and  makes  the  refinement 


sufficiently  strong. 
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6.  CONVERGENCE  OF  FEEDBACK  FINITE  ELEMENT  METHODS 

We  shall  refer  to  a  transition  operator  (and  the  corresponding  feedback 
f.e.m.)  as  5-regular  if  it  bisects  a_t  least  one  element  I  with  n  1(I’t) 

>  f  max  r;  *(I).  Since  the  transition  operators  A^^  through  A^^  all 

If  * 

i  A.  A. 

bisect  an  interval  I*  with  m  (I*)  =  max  "  (I)  these  are  indeed 

If  A. 

1-regular. 

Theorem  6.1.  Let  u^,  i  =  1,2...  be  the  finite  element  solutions  computed 
with  a  5-regular  feedback  f.e.m.,  0  <  6  S  1.  Then 

l-eillE  =  ^lu  ”  Ui^E  ’”  °  aS  1  ■*  °° 

provided  only  (2.2a-b)  and  (2.3)  hold. 

The  following  twTo  lemmas  are  essential  to  the  proof  of  Theorem  6.1. 

Lemma  6.1.  Let  H  be  a  Hilbert  space  and  let  -  S?  -  •  • •  ^  S^  -  -  • • •  _  H 

be  a  nested  sequence  of  closed  subspaces.  If  denotes  the  orthogonal  pro¬ 
jection  onto  S^  and  P  denotes  the  orthogonal  projection  onto 

CO 

S  =  U  S.  .  then  for  any  u  €  H 
i-1  1 

P . u  Pu  as  i  00 . 

l 

OO  ± 

Proof:  Since  Pu  f  U  S.  there  exist  v,  f  U  S.  =  S.  such  that 
-  ,  ,  i  i  .  ,  j  i 

i*l  j=l 


v^  -*  Pu  as  i 


-*>  O' 
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P.Pu  is  bv  definition  the  closest  element  to  Pu  in  S.,  and  thus 
1  1 

(6.1)  PiPu  -  Pu  as 

It  is  not  difficult  to  see  that  P.Pu  =  P.u,  and  this  in  combination  with 

i  i 

(6.1)  verifies  the  lemma.  0 


Lemma  6.2.  Let  H  be  a  Hilbert  space  and  let  H  L  2  S0  £  •••  IS,  I  ^i+1 *  ' ’ ' 
be  a  nested  sequence  of  closed  subspaces.  If  P^  denotes  the  orthogonal  projec- 

cx> 

tion  onto  S.  and  P  denotes  the  orthogonal  projection  ontoS  =  f|  S.,  then 
1  i=l  1 


for  any  u  (  H 


1 


t 

i 


P^u  ->•  Pu  as  i  -*  «. 

Proof :  Consider  the  sequence  s|  £  si,  £  • • *  £  £  •••  i  H.  The  ortho¬ 

gonal  projection  onto  s|  is  Id  -  P .  and  the  orthogonal  projection  onto 

— T  x  1  1 

U  S.  =  S  is  Id  -  P  (Id  denotes  the  identitv  operator).  The  previous 
i=l  1 

lemma  proves  that 


(Id-P^)u  -*•  (Id  -  P)u  as  i  -*• 

and  thus 


Piu  -*•  Pu  as  i  -*■  □ 

Proof  of  Theorem  6.1:  The  spaces  SP(A^)  £  S1^^)  -  •••  form  a  nested 
sequence  of  closed  subspaces.  The  elliptic  projection  P(A  )  is  the 
orthogonal  projection  onto  S' (i  )  in  the  space  H(7),  equipped  with 
the  innerproduct  B(-,*).  From  Lemma  6 . 1  we  conclude  that 


P(A1)  -  Pu, 
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OC 

P 

where  P  is  the  elliptic  projection  onto  'J  S  (A.),  i.e., 

i=l  1 


(6.2)  e.  =  u-u.  =  u  -  Pu  +  r . ,  :  0  as  i  -*■  «>. 

l  i  i  l 


As  before  let  R(A)  denote  the  elliptic  projection  onto  the  space 


{v  ^  H 1  ( I )  :  v(xp  =  0,  j  =  1 , .  . . ,N (a)  }  • 

According  to  (6.2)  and  Lemma  6.2 

(6.3)  R(Ai)e.  -*•  R° (Id  -  P)u, 


where  R  is  the  elliptic  projection  onto 


;  v  f 


H1  (7 ) 


v(x) 


0  at  all  nodal  points  x  of  <J  A.}. 

1=1  1 


For  each  i  >  1,  let  I.,..  (A.  be  an  element  which  is  bisected  in  the 

"  ’  j (i)  i 

or 

transition  from  A.  to  A. .  Since  U  A.  only  contains  finitely  manv  dif- 
i  i+1  .  .  i 

i=l 

ferent  elements  above  a  given  positive  size  and  since  repetions  in  the  sequence 

{1.,..}*  ,  are  excluded,  it  follows  that  h(I.,.J  -  0  as  i  -*•  °».  The 
j(i)  i=l  j(D 

corresponding  indicators 


‘<W 


,iR(7NIj(i))ei!E 


LJ I 


j(i) 


ll/2 


b(pi)' 


)dx 


with 
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.  =  R(A,)e. 

1  ii 


must  converge  to  0,  due  to  the  fact  that  h(l^^)  -*•  0  and  R(A^)e^  -*• 
R°(Id-P)u  in  H^(I).  The  feedback  f.e.m.  is  assumed  to  be  5-regular, 
0  <  5  5  1,  and  the  above  argument  may  thus  be  used  to  conclude  that 


(6.4) 


max  n  X(I)  -*•  0  as  i 


I«£. 

i 


From  (6.3)  and  (6.4)  we  get  that 


(6.5) 


R°  (Id  -  P) u  =  0. 


If,  to  the  contrary,  there  existed  a  point  x*  such  that  Ro(ld-P)u(x*) 

^  0,  then  | R(Ai)ei(x*) |  2  >  0  for  i  sufficiently  large,  because 

Ai  i 

of  the  convergence  (6.3);  let  (  A  denote  the  interval  that  contains 


x*,  then 


A  A 

>0  <  |R(Ai)ei(x*) |  <  Ch^/2  n  1(Iji), 


and  since  h.  =  h.  <1  this  would  implv  that 
J  J  - 


0  <  C  -  max  n  ( I ) , 

IfA . 


a  clear  contradiction  to  (6.4).  The  identity  (6.5)  in  combination  with 
(6.3")  and  (2.10)  leads  to 

t  (£  )  -*  0  as  i  -»  <*-, 
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and  an  application  of  Lemma  2.1  now  gives  the  desired  result.  □ 

Remark  6.1.  Note  that  Theorem  6.1  is  slightly  different  from  most  finite 

element  convergence  results  in  that  it  guarantees  convergence  without  any 

& . 

assumption  that  the  maximal  element  size,  h  ,  approaches  0.  Indeed 

if  u  is  a  piecewise  polynomial  of  degree  Sp  in  part  of  the  interval  I 
A . 

then  h  will  not  tend  to  0. 

Remark  6.2.  A  convergence  result  similar  to  Theorem  6.1  may  also  be 
proven  in  two  dimensions,  see  [  2  ]. 


2] 

7a.  QUASI-OPTIMALITY,  PRELIMINARY  RESULTS 

As  alreadv  pointed  out  in  section  3  we  distinguish  between  a  feedback 
f.e.m.  and  an  adaptive  feedback  f.e.m.;  an  adaptive  method  is  a  feedback 
method  which  is  optimal  with  respect  to  some  given  performance  measure. 

If  the  performance  measure  is  1  whenever  the  finite 
element  solutions  converge,  i . e . ,  whenever  u^  —  u  in  H  (7),  and  0 
otherwise,  then  Theorem  6.1  verifies  that  any  f-regular  method, 

o  ^ 

0  <  5  2  1,  is  adaptive  with  respect  to  this  measure  for  any  u  t  H  (I). 

In  this  section  and  the  next  we  shall  study  adaptivity  with  respect 
to  a  stronger  measure,  which  is  1  whenever  the  finite  element 
solutions  are  quasi-optimal  and  0  otherwise.  The  optimal  approximation 
is  defined  through  the  expression 

(7.1)  i(u,N)  =  inf  min  :u  -  v!£  , 

A:binary  mesh  -p,,. 

N(.:)=N  V€S  U) 

and  the  aforementioned  measure  will  be  1  iff 

(7.2)  liu  -  u  I!e  <  C<  (u,N(Ai))  , 

C  may  depend  on  u,  but  is  independent  of  i. 

Consider  the  following  simple  example: 

Example  7,1 

Let  u  i  H* ( 7)  \  { 0}  be  a  piecewise  polynomial  of  degree  SP  on  the 

mesh  0  =  xQ  <  xx  <  x2  =  1,  where  x.^  is  not  of  the  form  k/2J.  Then 

-  u  ^  0  for  anv  binarv  mesh,  but  at  the  same  time  i  (u,N)  =  0,  f°r  an>' 

l  E  " 

N  2  2,  and  hence  (7.2)  can  never  be  satisfied. 


i 


U 


As  is  evident  from  this  example  we  cannot  expect  to  verify  (7.2)  for 


any  of  our  feedback  f.e.m.  unless  we  require  that  1 (u,S)  not  become  too 
small.  We  shall  make  the  following  assumption  about  u, 


(A.  1) 


3  :  >  0  and  >  0  such  that 
J  (u,N)  2  C^N~  V  N  ?  1. 


The  indactor  r.  corresponds  to  the  finite  element  solution  of  (2.1a-b) 
and  it  is  therefore  not  entirely  determined  only  from  the  behaviour  of  u  in  IT. 
Assumptions  about  r)7  would  in  practice  be  very  difficult  to  verify  .  Instead  we 

A 

introduce  an  indicator  £7  which  corresponds  to  an  interpolant  of  u  and  which 

therefore  is  entirely  determined  based  on  knowledge  of  u  in  I^T;  in  the 

appendix  it  is  shown  that  the  ' ’ s  and  the  n's  are  intimately  connected. 

Let  u*(A)  =  Pa(£)u  €  SP(A)  denote  the  solution  to  (2.6)  in  case  a  =  1, 

d  A 

b  =  0.  Note  that  u^  on  the  interval  1^  is  simply  the  L^  projec¬ 
tion  of  u  onto  polynomials  of  degree  2  p  -  1  and  that  u*  Interpolates 
u  at  the  nodal  points  of  A;  u^  is  thus  locally  determined  from  u. 

Define 


(7.3) 


r  A 

"  j 


C(IA) 

1 


(^  (u-  u*(x)))2  dx 


i; 


=  mm 


q  LIA 


.  ,  (s»-’v 


dx 


1/2 


where  the  minimum  is  taken  over  all  polynomials  of  degree  5p  -  1. 


A  common  feature  of  the  transition  operators  A^  -  for 
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0  -  '  1  is  that  they  bisect  all  intervals  with  indicators  that  are  com¬ 

parable  (within  a  constant)  to  the  maximal  indicator.  One  might  expect 
that  this  will  lead  to  an  equilibration  of  the  set  of  indicators,  and  the 
intuition  is  now  that  this  equilibration  of  error-indicators  ensures  quasi¬ 
optimality.  In  order  to  make  this  plausibility  argument  rigorous  we 
naturally  have  to  make  the  assumption  that  there  exists  a  sequence  of  meshes 
for  which  a  certain  approximation  converges  at  optimal  rate  and  for  which 
at  the  same  time  the  errors  element  by  element  are  equilibrated.  Our  exact 
assumption  about  u  is 


(A. 2) 


For  any  N  .•  1  there  exists 
(i)  N(2(N))  =  N 


(N) 


such  that 


(ii) 


UI€. 


I  ( t (I) ) 2 
(N) 


1/2 


<  C . 1 ( u , N )  and 


(iii)  max  1(1)  2  C„  min  £(I) 


fc-  some  subset  of  elements 


7  (N)  .  (N) 


with  the 


(N)  ~(N) 

property  that  A  V  A  consists  of  at  most  K 

elements . 


The  constants  C ^  and  K  are  independent  of  N. 


Consider  the  following  example: 


Example  7.2 

Let  v^  :  I  ■*  1R  denote  the  function  given  by 
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V50  = 


for  0 


(-l)j  for  (j-i)2~k 


-k-1 


for  1-2 


-k-1 


(j  +^)2"\ 


1  £  j  5  2"  -  1 


rX 


u(x)  =  )  a,  '  v  (t)dt 

i0  K 


j  v  (t)dt 
J0 

fX 

!  V  (t)dt 

Jo 

rx 

v3(t)dt 


k=l 


square  summable. 


rX  a'  , 

7  a.  V  (t)dt  €  H 1(7), 
•'0  k-1 


2’ 

Let  4^  denote  the  mesh  ij*2  anc* 


let  u*  denote  u.(A.)  in  the  case  p  =  1.  It  is  not  difficult  to  see  that 

*  *  c 

the  v  ’s,  1  <  k  £  £  -  1,  are  constant  on  the  elements  of  £  and  that  the 

K  X 

v,  's,  c  S  k,  on  each  element  of  A  are  orthogonal  (in  L~)  to  the 

k  x 

constants;  we  thus  get 


d  .  1  x  r 

^(u-u*)  =  kltVk* 


(7.4) 


25 


The  right  hand  side  of  (7.4)  is  symmetric  in  the  nodal  points  of  L.  as 
a  consequence 


(7.5) 


/c 

j 


independently  of  j . 


From  (7.4)  and  the  fact  that  the  v^'s  are  indeed  orthonormal  in  L7(7), 
it  follows  that 


(7.6) 


J0 


kdx 


f  0 

(«-■£»- 


dx 


Pick 


ak  =  2_k  for  k  =  2j,  j  -  0,1,2,... 

a.  =  0  otherwise, 

k 


(7.6)  then  gives  that 


(7.7) 


f1  ,  d  , 

Jo  <E<U'“ 


*))2 


dx 


CN(A  )' 


£  =  2 


m 


For  C  =  2™,  we  define  £  '  =4  +  1  and  let  A1  denote  the  mesh  with  nodal 

^  £ 

points 


{<j-j>2  lK=1  U  {j2  }* 

for  which  clearly  N(aJ)/N(a^)  -*-1  as  m  -*  <*. 

exactly  so  that  the  v^'s,  k  =  2^  ,  0  5  ]  s  ro, 

and  from  the  selection  of  the  a,  ' 

k 


£  ' 

* 

=0 

The  mesh  A ^  is  constructed 
are  constant  in  each  ele- 
s  we  thus  get 


ment  of  A  ' , 
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We  may  add  N  ( 1 ' )  -  N(y  )  nodal  poincs  to  the  first  element  of  ^  with- 

'X  X 

out  upsetting  (7.7);  instead  of  (7.5)  we  now  get 


(7.9)  Cj5"  =  r  f'  for  at  least  values  of  j, 

with  v  arbitrarily  close  to  1  (provided  m  is  sufficiently  large).  It 

t 

follows  directly  from  (7.7)  and  (7.8)  that  the  u^  are  not  quasioptimal  even 
though  the  indicators  are  equilibrated  in  the  sense  of  (7.9).  c 


The  problem  with  the  preceding  example  is  not  that  i  dies  off  too 
fast,  rather  it  has  to  do  with  the  fact  that  there  is  a  total  lack  of 
correlation  between  the  length  of  an  element  and  the  associated  inter¬ 
polation  error.  To  avoid  this  problem  we  shall  assume  that 


(A. 3) 


■ 

There  exists  a  finite  set  of  point  { z ^ ^  and  a  constant 
C.j  such  that  V  I  =  (x,y)  £  7  : 

1 "  ■  * 

\ 

an  5  c3  min{ca1/2)  ,ca1/2)> 

where  I~^2  =  (x,^~^L)  and  X^2  =  (^^,  y)  . 


The  assumptions  (A.1)-(A.3)  are  easily  seen  to  be  satisfied  for  suf¬ 
ficiently  smooth  u  and  they  also  hold  for  a  wide  class  of  singular  u, 
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such  as  u(x)  =  xQ  -  x.  In  the  following  example  we  demonstrate  this  for 

Q 

the  case  u(x)  =  x  -  x  and  p  =  1. 


Example  7.3 

Let  u(x)  =  xa  -  x  €  ft3(7),  1/2  <  a.  It  is  not  difficult  to  see  that 


(A.l)  holds  with  c  =  1. 

Let  tp(x)  denote  the  function  cp(x)  =  x^2a  D/3  an(j  let  be 

the  mesh  with  nodal  points  x.  =  3(j/N)  =  (j/M)3/^2a  It  is  not  dif- 

2 

ficult  to  see  that  (£(!))  ,  1  <  j  <  K,  is  of  the  order 


2  h3 
J 


which  is  bounded  from  above  and  below  by  a  constant  times  N  .  Since 
i>(u,N)  2  C^N  1  this  leads  us  to  conclude  that  assumption  (A.  2)  is  also 
satisfied. 

Let  I  =  (x,y)  Q  1  be  any  interval;  it  is  not  difficult  to  see  that 
2 

(£;(!))  is  of  the  order 


K-^)2  u(y) I2  (y~K)3 

and  hence 

U(I))2  s  C|(^)2  U(y)|2  (y-x)3 

=  8C  min(!(^)2  u(^)|2,  |  (^) 2  u(y)|2}(^)3 
5  C’  min{f:(I1/,2)2, 


which  verifies  (A. 3). 


□ 


Definition  7.1.  Let  &  ,  i  =  1,2,...  be  a  sequence  of  meshes.  We  shall 

^  °° 

call  semi-equilibrated  if  there  exist  constants  C,  v  >  0  and 

subsets  c  a^  with  at  least  vN(ii)  elements,  such  that 

max  £,(1)  5  C  min  A  (I). 

Uh.  lit. 

1  1 


The  following  lemma  justifies  the  intuitive  feeling  that  equilibration 
implies  quasi-optimalitv. 

°1 

Lemma  7.1.  Assume  that  u  i  H  (7)  satisfies  (A.1)-(A.3).  If  a., 
i  =  1,2,...  is  a  semi-equilibrated  sequence  of  meshes  and  u^  =  u(A..)  is 
the  associated  finite  element  projection,  defined  in  (2.5)-(2.6),  then 

!  ei"E  =  '  U  ~  Ui‘‘E  “  CC(u,N(Ai)). 

Proof :  For  any  mesh  A  we  have 

u  -  u(A)l'p  sc  l  (e(i))2 , 

I€A 


where  C  depends  only  on  of  (2.2a-b).  It  thus  suffices  to  verify 

that 


l  (C(I))2  S  c  l  (c(D)2 

HA±  IEa’ 


where  A^  •  A  is  the  mesh,  of  same  cardinality  as  A^,  that  occurs 

in  the  assumption  (A. 2). 


Let  v,  A^  be  as  in  definition  7.1  and  let 


£  A^  consist  of 

those  elements  that  do  not  contain  any  of  the  points 


i 
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assumption  (A. 3).  /L  must  contain  at  least  vN(Ai)  -  K  elements  and 

since  it  is  only  necessary  to  prove  this  lemma  for  N(A^)  sufficiently 

large  we  may  assume  that  contains  at  least  v*N(A^)  elements, 

=  1c 

0  <  v'  <  v.  Divide  each  interval  of  into  2  subintervals  of  equal 

k 

size,  where  k  is  picked  such  that  2  v'>  1.  Since  there  are  at  least 
2^v'N (Ai)  2  N(A.)  of  these  subintervals  and  exactly  NtA.p  elements  in 
A ^ ,  it  follows  that  at  least  one  of  the  subintervals  Ip  is  entirely 
contained  in  an  element  Ip  €  A^.  From  the  definition  of  A^,  it  follows 


(7.10) 


max  £(I)  <  CC(J) 


-  CC3S(I0) 


s  cc^(l^). 


where  J  is  the  element  of  A^,  which  was  bisected  k  times  to  give 
Ip,  and  is  the  same  constant  as  in  (A. 3).  Based  on  (7.10)  and  the 

assumption  (A. 2)  we  conclude  that 


(7.11) 


C  (I)  £  CC(I') 


V  I  €  A.,  I*  6  A’ 


where  A^  =  A^^""i^  is  the  same  as  in  (A.  2).  Since  A^  and  A^  have 
a  comparable  number  of  elements  (7.11)  leads  to 


1  (cd))' 


C  J  (£(!))' 


S  C  y  (an)  , 


exactlv  as  desired. 
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7b.  QUASI-OPTIMALITY 

The  indicators  £  used  in  the  statements  of  the  preliminary  result 
in  section  7a  were  all  associated  to  the  interpolant  u^  (cf.(7.3)).  The 
transition  operators  A^^-A^^  defined  in  sections  4  and  5  were  all  based 
on  the  "actually  computed"  indicators  (cf.  (2.9)).  In  order  to  establish 
quasi-optimality  of  any  of  these  feedback  f.e.m.  we  thus  need  a  lemma 
relating  the  £, ' s  and  the  n's. 


Lemma  7.2.  Assume  that  a 
continuous  with  exponent 
such  that  for  any  mesh  t 


.  A. 1/2 

(V 


in  addition  to  satisfying 
1/2.  Then  there  exists  h 


0 


with  h^  S  h^  one  has 


C(hA)1/2  max  £(I) 
3  I€A 


nA  <  (aA)1/2  CA  +  C(hV/2 

1  III 


(2.2a-b)  is  Httlder 
(independent  of  u) 


max  i  ( I )  , 

UL 


t  i_l  -i 

1  5  j  S  N(a),  where  a.  =  a(-ji — r - J-)  and  C  is  independent 

J  - 

of  u  and  A. 


The  proof  of  this  lemma  is  given  in  the  appendix;  at  this  point  we 
proceed  to  verify  quasi-optimality  of  the  feedback  method  associated  with 
A^. 


Theorem  7.1.  Assume  that  a  is  Holder  continuous  with  exponent  1/2 

O  1 

and  satisfies  (2.2a-b).  Assume  that  u  €  H  (I)  satisfies  (A.1)-(A.3) 
d  v 

and  that  -r-  u  €  L  (I)  for  some  r  >  2.  Let  A.,  u . ,  i  =  1,2,...  be 
dx  l  i 

the  meshes  and  the  finite  element  solutions  generated  with  the  feedback 
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finite  element  method  based  on  the  transition  operator  ,  0  <  £  5  1. 

Assume  that  h  ^  S  h^  where  h^  is  as  given  in  Lemma  7.2.  Then 


ei'!E 


iu  -  u.!_  5  C4(u,N(A.)). 

1  L-  X 


Proof:  Let  I 


(1) 


denote  those  elements  of  L  .  which  are  bisected  in 

l 


passing  from  to  A^+^.  Since  there  are  at  most  finitely  many  different  ele¬ 


ments  in  U  above  a  given  positive  size  and  since  the  sets 

mutually  disjoint,  it  follows  easily  that 


(1) 


are 


(7.12) 


max  h(I) 
1 


as  i  • 


Let  I*  6  A .  denote  an  element  with 
l  l 


C(I£) 


=  max  t  ( I ) 

I€A. 

i 


Due  to  Theorem  6.1  we  know  that 


(7.13) 


£(!*)  0  as  i  -*•  <*, 


and  the  assumption  (A.l)  guarantees  that 


(7.14)  £(I*)  *0  Vi. 

A  combination  of  (7.13)  and  (7.14)  gives 

(7.15)  h(I*)  ->-0  as  i  ■*  ®. 
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Based  on  Lemma  7.2  and  (7.15)  we  conclude  that 


(7.16) 


i 

max  1(1)  5  C  max  n  (I) 

16*.  I€i. 

l  l 


for  i  sufficiently  large.  Similarly  Lemma  7.2,  (7.12),  (7.16)  and  the 


definition  of 


(1)  u  u  T‘i  ,  T(l) 

i .  show  that  for  anv  I.  cl. 
i  J  i 


(7.17) 


6'  max  1(1)  <  1(1 . 1)  , 

I€A.  3 

l 


with  0  <  p ' ,  provided  i  is  sufficiently  large. 

The  meshes  A  ,  i  >  ig  are  thus  according  to  (7.17)  generated  by 

^i 

always  bisecting  elements  that  satisfy  ?’  maxi(I)  5  £(I.  ).  Since  the 

UL.  3 

l 

indicator  1(1)  does  not  increase  through  bisection  of  an  interval,  we 
may  without  loss  of  generality  assume  that  each  transition  from  A  to 


"i-KL 


consists  in  the  bisection  of  only  one  interval,  satisfying  (7.17) 


(if  this  is  not  the  case  we  add  the  intermediate  meshes  to  the  sequence 
{A.}.  ,  ). 


We  divide  the  elements  of  into  three  distinct  categories 


those  I  £  A.  for  which  £(I)  <  8'C 


-1 


max 

IfAi 


1(1), 


i 


,  r 
"i 


those  I  €  A .  \  A® 
l  i 

point  {zk)£=1  of 


that  contain  one  (or  more)  of  the 
the  assumption  A. 3, 
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the  constant  A*  is  the  same  as  in  (7.17)  and  is  the  constant  of 

the  assumption  (A.  3).  We  shall  prove  that  contains  at  least  vNCa^) 

elements,  for  some  v  ^  0 .  According  to  the  definition  7.1  this  ensures 
that  the  sequence  of  meshes  is  semi-equilibrated  and  the  con¬ 

clusion  of  this  theorem  then  follows  directly  from  Lemma  7.1. 

Since  1  2  and  since  only  elements  that  satisfy  (7.17)  are  bisected, 

it  is  clear  that  the  element  bisected  in  passing  from  A.  to  ■' .  ,  is 

l  l+l 

S  IT 

not  in  A..  If  the  interval  being  bisected  is  in  L.  then  it  follows 
l  l 

from  (7.17),  the  assumption  (A. 3)  and  the  definition  of  aT  that  the  two 
resulting  new  elements  are  in  ,  i.e.,  we  obtain 

(7.18)  #Aj+1  5  Hl[  +  1 

{if  denotes  the  number  of  elements)  .  For  any  i^  +  1  <  i  define 

=  number  of  indices  j,  ifl  <  j  f  i-  1,  for  which  the  tran¬ 
sition  from  A.  to  A..,  consisted  in  bisecting  an 
J  J+1 

r 

element  of  A.. 

1 

If  we  prove 

M. 

(7.19)  lim  inf  —  >  v  •>  0, 

i-KXi 

then  we  immediately  get  from  (7.18)  that 

>  vN(A^),  for  i  sufficiently  large 

for  some  v  >  0,  which  completes  the  proof  of  this  theorem.  The  proof 


of  (7.19)  proceeds  by  contradiction:  assume  that 


M. 

— r — *•  0  for  some  (sub)sequence  i  -*■ 

this  implies  that  there  exist  arbitrarily  large  i  such  that  the  transition  from 

1.  to  -j+i'  ig  -  i  -  i  “  1  consist  at  least  i/2  times  in  bisecting  an 
2 

element  of  A..  Since  there  are  only  K  points  z^,  1  5  k  <  K,  in 

the  assumption  (A. 3)  we  conclude  that 


2  — T 

A.  contains  an  element  I!  of  size  2  ,  with  t  =  i/2K. 

l  l 


On  this  interval  I| 


«r.)2  <  !  A-'2 


,  (T*>  dx 


f  du  r  \r  ^~r 

<  i  i  (— )r  dx)  h ( i : ) 

j »  Q2C  /  1 


£  C2~Ul,  with  u  =  (l-|-)/2K. 


From  this  and  the  fact  that  I!  €  A.  it  follows  that 

l  l 


max  £(I) 2  £  C2-lJi  ,  or 


-uN(A.)  1/2 

i(u,N(A.))  £  C[N(Ai)2  1  ] 


for  certain  arbitrarily  large  i,  and  this  clearly  contradicts  the 
assumption  (A.l).  □ 
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Remark  7.1.  The  assumption  h  ^  can  obviously  be  replacad  by  the 

(weaker)  assumption  that  h  1  i  for  i  sufficiently  large.  This  is 

satisfied,  for  example,  when  the  solution  u  is  not  a  polynomial  of  decree 
ip  in  any  subinterval  of  I. 

Remark  1.2.  If  we  had  assumed  that  the  set  of  points  '  j<-l  was  emPty> 
i.e.  , 

C(I)  <  C3  min{£  ,  ■  (ljy2)  ' 

for  any  I  i  I,  then  the  set  L.  would  also  have  been  empty  and  we  had 
not  needed  to  assume  that  u  £  Lr(I),  r  >  2,  or  i(u,N)  >  CN  0  (only 

that  i (u,H)  "  0) .  □ 

Remark  7.3.  If  we  assume  that  the  set  of  point  of  the  assump¬ 

tion  (A. 3)  is  empty,  then  it  is  not  difficult  to  verify,  based  on  Lemma 
7.2,  that 


(7.20) 


A  . 

max  n  1(I)  s 

Ifi. 

i 


Ai 

C  max  h  (I) , 
ItV 


i  > 


i 


O’ 


where  o  1(I)  is  the  predictor  defined  in  section  5,  and  i  =  1,2,... 

(3) 

is  a  sequence  of  meshes  generated  by  the  transition  operator  A 
Consequently  there  exists  0  <  1  such  that 


c  * 


max 

^i 


n  1(I) 


), 


i  > 


i 


O’ 


for  any 


and  we  may  now  use  the  argument  of  the  previous  proof 


r 


to  show  that  the  finite  element  solutions  computed,  using  A 


,  are 


quas i-opt imal .  c 


(3) 


In  order  to  verify  quasi-optimality  of  the  feedback  f.e.m.  associated 

(2) 

with  the  transition  operator  A  we  shall  need  a  slightly  more  restrictiv 
assumption  than  (A. 3). 


(A. 3) ' 


There  exist  1  <  and  f  3  <_  1  such  that  for  any 

1  =  ( x ,  y )  c  1 


£(I)  £  C3  rain{:(I1^2),  r  ( I j )  -'  and 


max{£(l"/2),  f.(I1/2)}  £  3f  (I) 


,  X  +  v\  ,  ,+  ,x  +  v 

with  I^  =  (x ,  2  )  and  I ^  =  (-y^.y) 


We  could  permit  exceptional  points  as  an  (A- 3),  but  for  clarity 


of  exposition  we  have  omitted  these  here. 


Theorem  7.2.  Assume  that  a  is  Holder  continuous  with  exponent  1/2 

and  satisfies  (2.2a-b).  Assume  that  u  £  $^(I)  satisfies  (A.1)-(A.2) 

and  (A.3)'.  Let  1.,  u.  be  the  meshes  and  the  finite  element  solutions 
1  1 

generated  by  the  feedback  finite  element  method  based  on  the  transition 
C21 

operator  A'  ,  0  -  %  -  1,  0  <  ;'  <  1.  There  exists  0  <  Yq  such  that 

e.'ip  =  l:u  -  u  ;  £  C:(u,N(A  )) 

It  It  1 

provided  0  £  ^  (both  and  C  may  depend  on  u)  . 
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Assumption  (A. 3)'  implies  that  u  cannot  be  a  polynomial  of  degree 

■_  p  in  any  subinterval  unless  u  is  indeed  a  polynomial  of  degree  <p 

in  the  entire  interval  I;  this  latter  possibility  is  excluded  by  assump- 

2 

tion  (A.l).  From  Theorem  6.1  we  know  that  ,  r(I)~  -►  0  as  i  -<■  &, 

ie. 

l 

i.e.,  in  particular  max  "(I)  ~  0  as  i  —  <x.  Since  there  are  only  f initely  m?nv 

I  c  •  . 

'  l 

CL 

different  elements  in  i.  above  a  given  positive  size,  it  follows  that 

i=l  1 

either  h  1=  max  h(I)  -»  0  as  i  —  °c,  or  ;'(I)  =  0  for  some  subinterval 


1  •_  I,  1^0;  the  last  is  impossible  due  to  the  fact  that  u  is  nowhere 
a  polynomial  of  degree  5p,  and  we  thus  conclude 

(7.21)  h  1  =  max  h(I)  -  0  as  i  -*  cr. 


Theorem  7.2  will  be  proven  by  help  of  the  next  two  lemmas. 


Lemma  7.3.  Let  C.  denote  a  constant 

-  4 


C.  ->  max(C  ,1/6). 


Let  ip  be  a  given  index.  If  0  <  vq  is  taken  sufficiently  small  then 
there  exists  an  index  1^5  i  for  which  the  cardinality  of 


Af  =  (If  Ai  :  ay2£(I)  >  C^3  max  a^2r(I)} 

If  A. 

i 


is  at  least  N(£  )/2  for  any  0  £  v  £  (a^  denotes  the  value  of  a 

at  the  midpoint  of  the  interval  I) . 


- 


J. 


38 


Proof:  From  (7.21)  and  Lemma  7.2  we  see  that 


(7.22) 


if  If  and 

--  1(i)  >  „  max  1(I)  or  a^2C(I)  >  p  max  aj^fd),  P-  >  0, 
ia.  if- 

i  i 


then 

i'  ^  ( I )  =  a^/2Ul)  (1  +  PJ 


where  p.  -*•  0  as  i  -  independently  of  I.  Let  0  <  5  <  1  be  the 

(2) 

constant  entering  in  the  definition  of  A  and  select 

max  (e,<.'^)  <  5  <  1, 


then  it  is  clear  from  the  assumption  (A. 3)',  (7.22)  and  the  definition 

(2) 

of  the  transition  operator  A  that 


(7.23) 


max 
If  L 

i+1 


1/2.  ,,, 

aT  t(I)  < 


1/2., Ts 

max  a  ~ ( I) 
If  A  .  1 

l 


(2) 

after  each  application  of  A  ,  for  i  sufficiently  large.  Take  a  large  fixed 
S  i^  so  that  p i^  5  i,  are  sufficiently  small  and  (7.23)  is  satisfied  for 
i  5  i  (to  simplify  notation  we  assume  i^  =  =  1)  and  define 


Cn  =  max  a^2C (I) /min  a^2£(I); 

U  IfA  1  '  I€AX 

3  3 

if  Cq  s  then  there  is  nothing  to  prove.  If  we  consider 

the  following  sets  S^,  i  ?  1 : 

S.  =  fl  €  A.  :  aJ:/2Ul)  >  min  a?T/2T;  (I)  > ; 

1  11  If^  L 
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of  course  starts  out  to  be  A^ •  Define 


J.  =  max  (I)  /  min  a^“r(I) 

1  U*.  '  l€A, 


then  if  follows  directly  from  (7.23)  that 


(7.24) 


Ji+1  -  6Ji 


Let  k  >  1  be  the  first  integer  so  that 


J  «c  r 
Jk+1  -  L4 


Notice  that  k  £  k  ,  where  k  is  the  first  integer  to  satisfy 


0  3 

f  Cn  <  C  . 
0  -  4 


Since  J.  >  C. ,  1  £  i  £  k,  it  follows  that  whenever  the  transition  from 

i  4 

^i  to  A^+^,  1  £  i  £  k,  consists  in  bisecting  all  intervals  that 

satisfy 


A.  A . 

n  X(I)  ?  S  max  n  1(I) 
l€A . 


then  the  newly  created  elements  will  lie  in  S^+^.  Whenever  A''  bisects 
a  fixed  fraction  y  of  the  elements,  the  newly  created  intervals  may  all 
fall  outside  In  summary  we  get  that  for  v  sufficiently  small 


7 
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(7.25) 


S,  contains  at  least 
k+1 

(1  -  Y  -  (1  +  y)>-  -  -  (1  +  Y)  \) - ~~r 

(1  +  v) 


elements,  and 

J, 


k+1 


C4 


If  0  <  y  is  selected  appropriately  and  0  5  y  <  then  we  immediately 

see  from  (7.25)  that 


contains  at  least  N(y^+^)/2  elements,  and  2  ! 


this  verifies  the  lemma.  □ 

Based  on  Lemma  7.3,  it  is  now  easy  to  prove 
Lemma  7.4.  Let  C,  denote  a  constant 

-  4 

C4  >  max(C^,l/6) . 

There  exist  integers  i  and  k^  such  that  if  ®<Tq  is  taken  sufficiently  small 

and  0  5  v  S  Y-  then  the  cardinalitv  of 
u 


1/2  -3  3/2 

A.  =  {I  €  A.  :  aJ'cU)  >  CJ  max  a^c(I)} 

J  J  I€A. 

J 

is  at  least  N(a^)/2.  for  arbitary  (  €  N  and  some  i  +  k^  <  j  5  i  +  k^Cf+l). 

Proof :  Since  the  proof  in  many  ways  is  very  similar  to  the  previous,  we 

shall  omit  some  of  the  details. 

Pick  0  <  Yq  sufficiently  small  and  pick  i  such  that  the  assertion 


j 


I 


of  Lemma  7.3  holds.  J.et  k  be  the  smallest  integer  so  that 


(7.26) 


1/2 

max  a^  1(1) 


C.  ^  max  a (I) 

4  K£.  1 

1 


Dur  to  (7.23)  it  is  clear  that  k  <  k^  where  k^  satisfies 


5  0  5  c-1 

4 


The  size  of  0  <  v.  is  furthermore  taken  so  small  that 


(7.27) 


±  _  Y  -  (l+v)Y  - 


(l+v)k_1Y  >  0, 


0  5  Y  5 


For  is  j  S  i  +  k  define 


1/7  -4  l/'7 

S.  =  {If  L.'.a'lW  >  C  max  a  /‘KI)}; 

2  2  I€A. 


using  (7.22),  (7.26),  (7.27),  the  definition  of  i  and  the  assumption 


(A. 3),  we  get  that  for  0  S  r  S  Yr 


(7.28) 


S.  has  cardinality  at  least  N(4.)/2, 
J  J 


i  S  j  S  i  +  k . 


The  definitions  of  and  combination  with  (7.26)  gives 


Si+k  “  Ai+k  ’ 


and  the  desired  result  follows  now  directly  from  (7.28),  in  the  case 
(.  *  1 ;  the  case  of  £  >  2  follows  by  induction.  o 


Proof  of  Theorem  7.2.  From  Lemma  7.4  we  immediatelv  conclude  that  there 


exist  a  C  such  that  for  any  j 


l 
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i.  =  {I  €  2,.  :  C(I)  2  C  1  max  5(1)  } 

2  2  UL 

j 

CK 

has  cardinality  at  least  N(A.)/2.  According  to  Definition  7.1,  {A  } 

J  J  J 

is  therefore  semi-equilibrated,  and  the  theorem  follows  from  Lemma  7.1. 

□ 


8.  A  SIMPLE  COST  ANALYSIS 

Let  us  assume  that  the  cost  of  computing  the  function  u(£_),  given 
the  mesh  f  ,  is  proportional  to  N(a_.)  ,  >.  >  0  (the  cost  of  computing 

the  indicators  ^ . 1  and  generating  the  mesh  is  considred  negligibl 

compared  to  this).  The  total  cost  of  the  feedback  f.e.m.,  at  arriving  at 
u(a^)  ,  may  therefore  be  estimated  by 

i 

(8.1)  l  N(£.)\ 

j  =  l  3 

In  the  case  of  the  transition  operator  we  have  for  any  y'  <  v 

(and  N(£.)  sufficiently  large)  that 

N(;J  ( 1  +  Y  '  )  s  N(Aj+1), 

and  using  the  expression  (8.1)  we  get  the  total  cost  estimated  by 

i-1  . 

I  U+r'fJA  n(a,)a 

3-0 

y'  s  ■  N(A.)\  o  <  y'  <  y 

(1  +v')X-l  1 

(2) 

The  transition  operator  A  with  0  <  y  2  1  thus  has  the  property 
that  the  total  cost  is  comparable  to  the  cost  of  computing  the  last 
finite  element  solution.  The  same  is  in  general  not  true  for  the  tran¬ 
sition  operators  A^  and  A^:  if  the  function  u  has  a  strong  sin¬ 
gularity,  then  a  few  indicators  near  this  singularity  will  remain  larger 
than  the  rest  through  successive  refinements,  and  the  total  cost  will  be 
of  a  quite  different  magnitude  than  the  cost  of  a  single  finite  element 
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solution.  In  the  program  FEARS  this  problem  has  been  solved  for  the  tran- 
(3) 

sition  operator  A  through  the  introduction  of  the  aforementioned  short 
passes  (cf.  section  5). 


/ 
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9.  APPENDIX 


°1 

The  bilinear  form  B  introduced  in  (2.4)  coerces  the  norm  on  H  (')> 

,1/2 


actually  bv  definition 


.u!  £  =  [B(u,u)  ]' 


Due  to  the  form  of  B  it  is 

°1 


also  very  natural  to  consider  it  relative  to  the  normed  spaces  W  (7 )  and 

W^(7).  As  the  norm  on  W^(7)  we  choose  j i u { ! _  =  sup!-^-  u(x) I  +  sup'u(x)t 

1  “  1.“  ..rVdx  1  1 


O 

and  as  the  norm  on  W^(7), 


1,1 


,1 


0 


x€l 

1 


—  u  I  dx  + 
dx  1 


u  dx. 


Theorem  9.1.  Let  a  and  b  satisfy  (2.2a-b).  Assume  that  a  is  Holder 
continuous  with  exponent 
that 


0  <  6  S  1.  There  exist  6  and  h^  such 


inf 

ufSP(A) 
lIU'!l,cc  =  1 


sup  jB(u,v) 
v<?Sp 


>  X 


0 


for  anv  mesh  with  fi  $  . 


Before  we  give  the  proof  of  this  theorem,  let  us  show  how  the  result 
implies  Lemma  7.2. 

Let  u^(A)  =  P*(A)u  €  SP(A)  denote  the  solution  to  (2.6)  in  the 
special  case  a  =  1,  b  =  0  (as  in  section  7).  If  e(A)  denotes 
u  -  u(A)  and  e^(A)  =  u  -  u^(A) ,  then  according  to  (2.9)  and  (7.3) 


nj  =  I!R(I  NiJ)e(A)|lE 


.A 

"j 


e.C))2  dx 


,  1/2 


For  any  v  -  SP(A)  we  have  that 


f 


B(u(a)-u*(a) ,v) 
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B(eA(4)  -  e(A),v)  =  B(e;f.(A),v) 


(9.1) 


N(A)  r  ,  ,  N(A) 

[  |  a  —  e  (A)  —  v  dx  +  / 

j-i  i,a  d*  *  dx  j-i  ‘j. 


bei(4)v  dx 


1A 

1 


N(A) 


I  j  ^"Vd^  e*(A)  d^vdx  +  /  j  be^  (A)v  dx 


j 

N(A)  (■ 


j  =  x  ^ 


j=l  J  jA 


A  A 

^  x;,+x.  . 

where  aT  =  a(— ^ ^ ;  in  (9.1)  we  have  used  the  fact  that  e  (A)q  dx 

J  1  J  a  dx  * 


A 

I . 
J 


0  for  any  j,  and  any  polynomial  q  of  degree  Sp  -  1. 


Since  a  is  Holder  continuous  with  exponent  1/2,  la  -  a^  *  5  *VbT)  on  * 
so  that  based  on  (9.1) 


|B(u(A)  -  u*(A),v) 


[  N(a) 


-c(  l  sup  1^  v(x)  I  I  |— -e^(A)|dx+  l  sup  |v(x)  I  j  |e*(A)|dx 

\j  =  1  vfIA  TA  ^  =  1  xfIA  ’  IA 


N(A) 


(9.2) 


'N(A) 

5  c(  I, 

.J  =  l 


N(A) 

sc  l 

j-1 


xfl. 


|fv|dx/f  -d 
A  dX  \  j  a 


IA 

J 


j=1 


(r  .  n  \  1/2  N(a)  a  ,  f 

J^e*(A)l  dxj  +  J1(hj)  J  lv,dx 


IA 

J 


j 


I  d  |  ,  ..A.  1/2 

—  v  dx  +  (h.) 
dx  j 


|e  (A)  ; dx 


A 

IT 

1 


A  i v  !  dxj 

h  / 


‘i- 


In  the  last  inequality  we  have  used  the  fact  that  e^(A)  "0  at  the  nodal 
points  of  A .  it  follows  now  directly  from  Theorem  9.1  and  the  estimate 


(9.2)  that 
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u(a)  -  u  (a)!' 

K  1  ,C 


C  max  £ (I) 
1 64 


Hence 


tk  =  jlR(lMA)e(A)l!r 
J  J  t. 


a(-^  (u-u^(a)))2  dx 


b  (u-u^(a) )  dx 


a^dx  (u*(a)-u(a)))2  dx  +  ^  b  (u.v(A)-u(A) )  2  dx 

if  l" 

J  J 


?  (aV/2^ 

J  J 


C(h'T) 1/2  max  £  (I) , 
J  l€A 


and  similarly 


ru  5  (a“)1/2  +  C(t01/2  max  £(I)  . 

2  2  2  2  IfA 


In  our  proof  of  Theorem  9.1  we  use  the  following  observations. 


Lemma  9.1.  Let  v  be  a  polynomial  of  degree  5  r.  Let  I  =  (y,z)  and 
a(x)  €  Lro,  0  <  <  a(x)  <  <  <* .  Assume  that  v  is  given  by  the 

expression 


l  c.(x-y)1 


Define  the  norms 
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-  2  2i 

;  c  h 

i=0 


Tl/2 


I  !ci!hl 

_  r\ 


i=0 


vi""  =  A 


r  |f  a(x) v(x) (x-y) idx 

r*  '  *  i. 

I  - 1 - 

i=0  h 


with  h  =  z-y.  There  exists  a  constant  C  independent  of  I,  v,  a, 
(but  depending  on  and  r)  such  that 

1/Cl'vl!"  s  |lv||L  (i)  s  C|!v|" 

cc 

and 


I!  vi 


L^I)  " 


cl:  vij 1 


The  proof  of  this  lemma  consists  of  a  simple  scaling  argument  and  the 
fact  that  all  norms  are  equivalent  on  a  finite  dimensional  space.  That 
C  depends  only  on  and  r  follows  directly  from  this  argument. 


Proof  of  Theorem  9.1.  Let  us  first  consider  the  case  that  b  =  0  in  the 
definition  (2.4)  of  B. 


Given  u  €  Sf(a)  let  I  =  (x.  x.  )  be  an  interval  with  the 

J0  J0 

property  that  the  maximum  of  uj  is  attained  in  I.  Let 


P_1  i 

L  c  (x-x  )  , 

i=0  J0  1 


x  €  I, 


and  define 
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P-1 

z(x)  V  d,(x-x.  )1, 

i=0  1 


x  f  I 


where  the  d/s  are  selected  so  that 


a (x) z(x) (x-x .  )i  =  h"i  c.. 

:  I  V1 


0  <  i  <  p  -  1 


Using  Lemma  9.1  we  thus  get 


(  a(x)  -j—  u(x)z(x)dx 


I 


dx 


(9.3) 


P-1 


)  c  a(x)z(x) (x-x.  )  dx 

i-1  1  J  I  V1 


P-1 

j.  c;h 


i=l 


2i  d  .  2 

-  fi:dJ  U  L  (I) 


dx  U''l  (T)  ’ 


and  furthermore 


,2|lLl(D  5  C!iZl:L’(I) 


(9.4) 


=  c  y 

i=0 


P-1 


,_1  j / j.a(x) z(x)  (x_xj  _j)ldxi 


c  l  KJh1  =  C|!^  u|! " 

i-0 


5  Cit-r—  ut!_  .  T . 
dx  L  (I) 


C  —  u !' 

dx  "L  (I) 


Extend  z  to  be  zero  outside  of  I  and  define 
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:  (x)  = 


fX  1 

A  1  dt 

’ o  a(t) 


where  A  is  chosen  such  that 


,1 

it  ( 1 )  =  z(x)dx 

Jo 


It  is  clear  that 


C  z 


L 1  ( I )  ‘ 


Let  v  be  the  piecewise  linear  iijrterpolant  of  c  on  £ ,  then 


(9‘5)  'dx  l  (I)  ~  c!Aj(h“>‘  -  C  z,  (j)  (h-) 


S  C!'^U::L  a)^9* 


Define 


v(x) 


z(t)dt  -  v (x) , 


it  is  clear  that  v  €  Sp(i)  (the  construction  of  v  guarantees  that 
v  =  0  at  the  endpoints  of  I),  and 


d  d  , 
a  —  u  —  v  dx 
dx  dx 


rld  f1  d  d 

a— uzdx-  |  a  —  u  -z—  c  dx 
0  dx  J0  dx  dx 


,] 


d  .  d  .  d  .  , 
a  -t~u(  —  v  -  — a;)dx 
dx  dx  dx 


0 
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ri  ,  -1  . 

d  .  id 

a  —  u  z  ax  -  A  —  u  dx 

■  o  dx  .fo  dx 


0 


d  ^  r 

3  d^d^  (v-:)dx. 


Due  to  the  estimates  (9.3),  (9.5)  and  the  fact  that  u(0)  =  u(l)  =  0,  we 
now  get 


(9.6) 


c1  d  d 

x  ,jL 

„2 

;  a  —  u  — —  v  > 
j q  dx  dx 

dx 

U  L  (I)  - 

(X 

> 

u!:2 

1  dx 

r  d  2 

C  dx  U  L  (I)(h  > 


t  ru! 


,2 


1,- 


fox  h1"  sufficiently  small.  From  (9.4)  and  (9.5)  we  immediately  see 


v !  <  p  i  — — 

1,1  dx  V  L1(7) 


C(  Z:L1(I)  +  dx  V'li(1) 


-  cilT"  /t^  5  c;u'i  , 

dx  L  ( I)  1,°" 


and  a  combination  of  this  and  (9.6)  leads  to  the  desired  result  in  case 
b  =  0. 

We  now  consider  the  general  case  of  a  nonzero  b.  Let  v  be  the 
same  function  as  before  and  set  =  v  +  w,  with  w  to  be  picked  later. 


Then  we  have 
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B(u,v1) 


f1  (J  d  f1 

,  a  —  u  —  v  dx  +  b  u  v  dx 

Jo  dx  dx  Jo 


rl  d  d  rl 

+  |  a  — u  — wdx  +  buw 

Jo 


dx  dx 


dx . 


Let  w  denote  the  solution  to  the  two  point  boundary  value  problem 


-r~(a  -7—  w)  +  bw  =  -bv 
dx  dx 


w(0)  =  w(l)  =  0  ; 


it  is  not  difficult  to  see  that 


|^w(x)  -  ~V(X+C)j  <  C-v!^  c9 


5  C  u: 

1,® 


We  pick  w  to  be  the  linear  interpolant  to  w  on  the  mesh  A,  then 
=  v  +  w  €  SP(A)  and 


B(u,v  )  >  [  a  v  ”  Cjiui!  |'w  -  w|’ 

1  Jq  dx  dx  1,®  l,c 


(9.7) 


2  5!lu|l2  -  Cllujl.2  (hA)6 

1,®  1,” 


for  hA 


sufficiently  small. 


Similarly  as  before  we  also  have  that 


(9.8) 

and  a  combination  of  (9.7)  and  (9.8)  leads  to  the  desired  result.  □ 


f 
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